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The assumption that our Universe is close to a late time fixed point of the equations of cosmol- 
ogy, leads to a modification of the latter to include energy exchange between the matter and the 
"dark energy". The brane- world scenario provides a natural set-up for such energy exchange and 
is analyzed in detail. The role of brane-bulk energy exchange and of an induced gravity term on 
a single braneworld of negative tension and vanishing effective cosmological constant is studied. It 
is shown that for the physically interesting cases of dust and radiation a unique global attractor 
which can realize our present universe (accelerating and Q<Q. m o < 1) exists for a wide range of the 
parameters of the model. For fl m o =0.3, independently of the other parameters, the model predicts 
that the equation of state for the dark energy today is wde.o = —1-4, while Q m o =0.03 leads to 
wdb, o = — 1.03. In addition, during its evolution, wde crosses the wdb = —1 line to smaller values. 



I. MOTIVATION - INTRODUCTION - 
CONCLUSIONS 

The use of renormalization group and fixed point ideas 
in various efforts to understand the symmetries or the pa- 
rameter values chosen by Nature, is not new 0- Imagine 
that the equations of cosmology have a late time fixed 
point. This means that as long as the Universe is old 
enough, the values of its observable quantities, simple ex- 
pressions in terms of the fundamental parameters, should 
be close to their values at the fixed point. Then, the an- 
swer to the coincidence question "why it is that today 
f2 m o and o e,o are of the same order of magnitude" , in- 
stead of relying on a fine-tuning of initial conditions, it 
reduces to an appropriate choice of the parameters of the 
theory, which should not require any serious fine-tuning. 

However, as we will now demonstrate, the equations of 
standard cosmology do not accommodate our Universe 
as a late-time fixed point. Specifically, I will show that 

If the energy density of a perfect fluid with equation of 
state w>— 1/3 of any cosmological system is conserved, 
all fixed points of the system with Q m ^ are decelerating. 
This contradicts the acceleration of our Universe. 

Indeed, with p the energy density of the perfect fluid, 
the equations of cosmology are 

H 2 -2 1 {p+p DE ) = 
p + 3(1 +w)Hp = 
Pde + 3(l + w DE )Hp DE = (1) 

where 7 = AttGn/3- The first is the standard Hubble 
equation. The second is the matter energy conserva- 
tion. The equation for p EE can always be brought into 
the above form, where wde is in general time as well 
as model dependent. Using QJ, it is straightforward 
to derive d(Cl m /QDE)/dhia — 3(fl m /flDE)(wDE — w) 



and 2q = l + 3(wfl m +WDE^DE), where fi m = 2jp/H 2 , 
Ude = ^1PdeI 'H 2 and q = —d/aH 2 . At the fixed point 
(denoted by *) d(n m /flDE)/dlna — 0. For Q m *^0 one 
obtains w EE * — w i and consequently, 2q« = l + 3w > 
Q.E.D. 

Thus, independently of the cosmological model, the 
only way our accelerating universe with f2 m * ^ can be 
close to a late time fixed point is by violating the standard 
conservation equation of matter. 

Braneworlds are a natural set-up for such matter non- 
conservation. The observable Universe is represented by 
a 3-brane embedded in a higher dimensional bulk and 
the above violation may be the result of energy exchange 
between the brane and the bulk. In particular in five di- 
mensions, a universe with fixed points characterized by 
O m * 7^ 0, q* < was realized in Q in the context of 
the Randall- Sundrum braneworld scenario with energy 
influx from the bulk. However, those fixed points cannot 
represent the present universe, since they have f2 m * > 2. 
Here I will present a brane-bulk energy exchange model 
with induced gravity, whose global attractor can repre- 
sent today's universe. Four-dimensional scenaria with 
accelerating late time cosmological phase may be found 
in 0. 

An interesting result: Let us, then, generalize the equa- 
tions of cosmology to allow for violation of matter energy 
conservation, due to some hypothetical "interaction". 

H 2 -2 1 {p+ PDE ) = 
p + 3(l + w)Hp = -T 
Pde + 3{l + w DE )Hp DE = T (2) 

Given the second equation of |J2J, the equation governing 
Pde can always be brought into the above form, where 
Wde is time and model dependent. If, in addition, at the 
fixed point 

ff*T*^0, p = pde = 0, (3) 
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one obtains 

wde* = -1 - 



1 + w 
0=1 _ 1 



(4) 
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A few comments are in order here, (a) Equation J2J is 
universal. It does not depend on the form of T or the 
function woEif)- (b) For f2 m * < 1 equation J3J gives 
wde* < — 1- (c) In the phenomenologically interesting 
case of w = and f2 m * = Qcdm = 0.3, i.e. if we take 
the present value of matter density to include the cold 
dark matter, one obtains wde* = —1-4. Finally, (d) 
assuming that the dark matter is not part of Pmatten 
n m * — ftbar — 0.03, and one obtains Wde* = —1-03. It 
is not surprising that once you allow for influx from the 
bulk, you may obtain constant acceleration in the ob- 
servable Universe. After all, this is similar to the steady 
state universe. What is surprising is that leads to the 
specific universal prediction for the dark energy equation 
of state, which in addition, has been favored by some 
recent analyses of the cosmological data. 

A braneworld model with a 3-brane embedded in a 5- 
dimensional bulk will be presented, whose equations will 
be explicitly brought into the form of equations J5J . For 
a wide range of parameters they have a late time global 
attractor, representing an accelerating universe with con- 
stant energy density on the brane, sustained by energy 
influx from the bulk. As we argued above, the fixed point 
satisfies ©, (0J. Moreover, via a numerical integration it 
will be shown that for a wide range of initial conditions 
the universe, during its evolution, crosses the wde = — 1 
line from higher values. Interestingly enough, these fea- 
tures are favored by several model-independent [j| as well 
as model-dependent 0, S IE 13 analyzes of the astronom- 
ical data. 
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where primes indicate derivatives with respect to y, and 
dots derivatives with respect to t. The five-dimensional 
Einstein equations take the usual form 



Gac — 



1 



tT 



AC\tot, 



2M 3 
where 

is the total energy-momentum tensor, 

diag(-A, -A, -A, -A, -A) 



Tq\v,B 

TX,b = diag(-V, -V, -V, -V, 0)^M 



Tc\ m ,b 



diag(-p,p,p,p,0) 
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II. THE MODEL 

We shall assume that we live on a 3-brane embedded in 
five dimensional AdS space-time. We consider the model 
described by the gravitational brane-bulk action [Iol| 



d 5 3 



^(M 3 R-A) 



^-h{m 2 R~V), (5) 



where R, R are the curvature scalars of the bulk metric 
gAB and the induced metric Hab — 9ab — nA n B respec- 
tively (n A is the unit vector normal to the brane and 
A,B = 0,1,2,3,5). The bulk cosmological constant is 
A/2M 3 < 0, the brane tension is V, and the induced- 
gravity crossover scale is r c = m 2 /M 3 . 

We assume the cosmological bulk ansatz 



ds 2 



-n(t, y) 2 dt 2 + a(t, y) 2 -f lJ dx t dx j + b(t, yfdy 2 , (6) 



where jij is a maximally symmetric 3-dimensional met- 
ric, parameterized by the spatial curvature k = —1,0, 1. 
The non-zero components of the five-dimensional Ein- 



TQ\m,B is any possible additional energy-momentum in 
the bulk, the brane matter content T^\ m ^ consists of a 
perfect fluid with energy density p and pressure p, while 
the contributions arising from the scalar curvature of the 
brane are given by 
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Assuming a Z2 symmetry around the brane, the singu- 
lar part of equations (|llf> gives the matching conditions 
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(the subscript o denotes the value on the brane), while 
from the 05, 55 components of equations l|ll|l we obtain 
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n p \i kb 2 a = T 55 -Ab 2 Q 

~nj\ a 2 6AP ' 
(21) 

where T05, T55 are the 05 and 55 components of T^c\m,B 
evaluated on the brane. Substituting the expressions 
(fT£|l . (fE?f in equations (J22J, J2J, we obtain the semi- 
conservation law and the Raychaudhuri equation 
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where iJ D = a /a n is the Hubble parameter of the 
brane. One can easily check that in the limit m — * 0, 
equation 1)23|) reduces to the corresponding second order 
equation of the model without R 3J. Energy exchange 
between the brane and the bulk has also been investi- 
gated in mmm. 

Since only the 55 component of TAc\m.B enters equa- 
tion (|23|l . one can derive a cosmological system that is 
largely independent of the bulk dynamics, if at the po- 
sition of the brane the contribution of this component 
relative to the bulk vacuum energy is much less impor- 
tant than the brane matter relative to the brane vacuum 
energy, or schematically 



Tl 



A 



(24) 



Then, for |A| not much larger than the Randall-Sundrum 
value V 2 /\2M 3 7 the term T| in equation (|23*|) can be 
ignored. Alternatively, the term Xg can be ignored in 
equation 123|) if simply 



If 



A 



< 1. 



(25) 



Note that relations (|24|l and (|25|l are only boundary con- 
ditions for T5 , which in a realistic description in terms of 
bulk fields will be translated into boundary conditions on 
these fields. In the special case where l|24l) . Q25JI are valid 
throughout the bulk, the latter remains unperturbed by 
the exchange of energy with the brane. 

One can now check that a first integral of equation 1|23|) 

is 



2H 2 (P +V | 6 
2m 2 r 2 
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with x satisfying 
X + 4n H x 



m 2 b n V 



p+V 
6m 2 



(27) 



and T = 2T® is the discontinuity across the brane of the 
05 component of the bulk energy-momentum tensor. The 
solution of iPffil for H a is 
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and equation (|27(l becomes 
x + 4n H x=^{l± ' 
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At this point we find it convenient to employ a coor- 
dinate frame in which b = n = l in the above equations. 
This can be achieved by using Gauss normal coordinates 
with b(t, z) = 1, and by going to the temporal gauge on 
the brane with n Q — l. It is also convenient to define the 
parameters 
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For a perfect fluid on the brane with equation of state p = 
wp the cosmology on the brane is described by equations 
PSty. (|2EJl, (EH, which simplify to (we omit the subscript 
o in the following) 



p + 3(1 + w)Hp = -T 



H 2 = p + 2 1P ± /3v/A + 24 7 p+x - -T 



4tf X - 24 7 T (l ± ±- VA + 24 7 p+ X ) , 



(34) 
(35) 
(36) 



while the second order equation Q23[) for the scale factor 
becomes 

a , . A + 6(1 — 3w) 7 p 

-=p-(l + 3whp±f3 2 JE. 37 
a v A + 24 7 p + x 



Equivalently, setting ip = ^/X + 24 7 p + x, equations 
JlU, ©, (123 take the form 



if 2 



4> + 2H 



p + 2jp±f3ip- 
A + 6(l-3w) 7 



k 



± 



2 7 T 



- = fi-(l+3w)jp±l3 
a 



A + 6(l-3w) 7 p 



(38) 
(39) 
(40) 



Throughout, we will assume T(p) = Ap v , with v > 0, A 
constant parameters 0, 0] . Notice that the system of 
equations l|34|) - (|36|) has the influx-outflow symmetry T — > 



4 



— T, H — > —H, t — > — t. For T = the system reduces to 
the cosmology studied in [l5| . 

We will be referring to the upper (lower) ± solution as 
Branch A (Branch B). We shall be interested in a model 
that reduces to the Randall-Sundrum vacuum in the ab- 
sence of matter, i.e. it has vanishing effective cosmo- 
logical constant. This is achieved for p, = =f0y/\, which, 
given that m 2 V+12M 6 is negative (positive) for branches 
A (B), is equivalent to the fine-tuning A= — V 2 /12M 3 . 
Notice that for Branch A, V is necessarily negative. Cos- 
mologies with negative brane tension in the induced grav- 
ity scenario have also been discussed in 0] . 

Consider the case k = 0. The system possesses the 
obvious fixed point (p*, i?*, tp*) — (0,0, \/X). However, 
for sgn(H)T < there are non-trivial fixed points, which 
are found by setting p — ip — in equations l|34|). 
For w < 1/3 these are: 

2T(P \ ] 2 2 =2/*+(l-3tt,) 7 p« 



9(l+w) 2 p 2 



± V '9(1 + w) 2 ^ 2 p 2 + 4/3 2 [A + 6(l-3iu)7p„] (41) 



(42) 



a ± 3(1 +w) + g _ 3u; ) ] = _ (43) 

P 

Equation lf4*f]|) gives 

'S\ T(p„) 2 



a/* 9(l+w) 2 p 2 ' 



(44) 



which is positive, and also, it has the same form (as a 
function of p*) as in the absence of R. The deceleration 
parameter is found to have the value 



q* 



-i. 



(45) 



which means i3* = 0. Furthermore, at this fixed point we 
find 



Equation (|41|) , when expressed in terms of VL m , 
one root for each branch 

6(l-3n;)/3±'/A(l-3w-4a-l; 
P * ~ 27 (2fi m l + l + 3w)(fi m l-l) 



(46) 



has only 



(47) 



However, it can be seen from (|47f) that for — 1 < to < 1/3 
and O m * < 1 the Branch B is inconsistent with equation 
(|4ip. On the contrary, Branch A with — 1 < w < 1/3 
and f2 TO * < 1 is consistent for < 6(1 — 3w;)/3+VA(l — 
3^-40-1) < 3V4(l-3w) 2 /? 2 -(l + u;) 2 A. Thus, since we 
are interested in realizing the present universe as a fixed 
point, Branch B should be rejected, and from now on 
we will only consider Branch A. So, we have seen until 
now that for negative brane tension, we can have a fixed 



point of our model with acceleration and < fl m * < 1- 
This behavior is qualitatively different from the one ob- 
tained in the context of the model presented in Q (for 
— 1/3 < w < 1/3), where for positive brane tension we 
have fi m * > 2, while for negative brane tension the uni- 
verse necessarily exhibited deceleration; therefore, in that 
model the idea that the present universe is close to a fixed 
point could not be realized. 



III. CRITICAL POINT ANALYSIS 

We shall restrict ourselves to the fiat case fc = 0. In or- 
der to study the dynamics of the system, it is convenient 
to use (dimensionless) flatness parameters such that the 
state space is compact 01 • Defining 



2 7 p 



H 
D 



where D = ^/H 2 — /i, we obtain the equations 



u m + City — 1 



(48) 
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l-UJr, 



(50) 



Z' = {l-Z 2 ) (1-Z 2 ) 



l-Z 2 -3(l-3w)[3 2 p- 1 uj. 



1—LJr. 



- 1 



l + 3u> 



,(51) 



with ' = d/dT = D~ 1 d/dt. Note that -1 < Z < 1, while 
both lo's satisfy < w < 1. The deceleration parameter 
is given by 



1 

~ZT< 



1 + iw 



(1-Z 2 



-3(l-3«j)/? 2 /i - V) 



1— w* 



(52) 

andi?' = -HZ(q+l). The system of equations lf5U|l - l(5T)) 
inherits from equations (|34|I - H36() the symmetry A — > —A, 
Z — > —Z, t — > —t. The system written in the new 
variables contains only three parameters. However, going 
back to the physical quantities H , p one will need specific 
values of two more parameters. 

It is obvious that the points with \Z\ = 1 have H = oo. 
Therefore, from (|38|l it arises that the infinite density p = 
oo big bang (big crunch) singularity, when it appears, is 
represented by one of the points with Z= 1 (Z = —l). The 
points with uj m = 1, \Z\^=1,0 have u>' m = oo, Z 1 = oo and 
finite p, H; for w < 1/3, one has in addition a/a= +oo, 
i.e. divergent 4D curvature scalar on the brane. 

The system possesses, generically, the fixed point 

(a) (uj m *, oty*, Z*) — (0,1,0), which corresponds to the 
fixed point (/?„, if*, = (0,0, VA) discussed above. 
For f < 3/2 there are in addition the fixed points 

(b) (w TO *,City*,.Z*) = (0,1,1) and (c) (ow, Ctty*, Z*) = 
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(0,1,-1). All these critical points are either non- 
hyperbolic, or their characteristic matrix is not defined 
at all; thus, their stability cannot be studied by first 
order perturbation analysis. In cases like these, one 
can find non-conventional behaviors (such as saddle- 
nodes and cusps 0|) of the flow-chart near the criti- 
cal points. There are two more candidate fixed points 
(d) (a; m *,w^*,Z*) = (1,0,1) and (e) (w TO *, w^*, Z*) = 
(1,0,-1), whose existence cannot be confirmed directly 
from the dynamical system, since they make equations 
|[5DJ, (|5"T)) undetermined. Apart from the above, there 
are other critical points given by 



2 7 



3(1 + w) Z* 



(53) 



6f3 2 



(1+3w)cj^+(1-3w) l-^-(l-Z 2 ) ^-2[l-(l-Z 2 ) 2 } 
L p J 

= 0. (54) 

They exist only for AZ* < and correspond to the ones 
given by equations (|41|) - (|43|l . For the physically interest- 
ing case w = with influx we scanned the parameter space 
and were convinced that for v ^ 3/2 there is always only 
one fixed point; for v<3j2 this is an attractor (A), while 
for v>3/2 this is a saddle (S). For u> = 0, 2/ = 3/2 there 
is either one fixed point (attractor) or no fixed points, 
depending on the parameters. For the other character- 
istic value w = 1/3, we concluded that for v < 3/2 there 
is only one fixed point (attractor) , for v > 2 there is only 
one fixed point (saddle), while for 3/2 < v < 2 there are 
either two fixed points (one attractor and one saddle) 
or no fixed points at all, depending on the parameters. 
For w = 1/3, v — 3/2 there is either one fixed point (at- 
tractor) or no fixed points. Finally, for w = 1/3, v = 2 
there is either one fixed point (saddle) or no fixed points. 
These results were obtained numerically for a wide range 
of parameters and are summarized in Tables 1 and 2. 





v < 3/2 


v = 3/2 


v > 3/2 


No. of F.P. 
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Nature 
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S 



Table 1: The fixed points for w=0, influx 
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No. of F.P. 
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A,S 
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Table 2: The fixed points for w=l/3, influx 

The approach to an attractor described by the linear ap- 
proximation of (|50|I - H51|I is exponential in r and takes 
infinite time r for the universe to reach it. Given that 
near this fixed point the relation between the cosmic time 
t and the time t is linear, we conclude that it also takes 
infinite cosmic time to reach the attractor. 

Defining e = sgn(H), we see from (|5T71l - <|5T|) that the 
lines Z = e (v < 3/2), ui m — are orbits of the sys- 
tem. Furthermore, the family of solutions with Z ss e 



and dZ/duj m — Z' /uj' m ss is approximately described for 
v < 3/2 by ui' m = e(l + 3w)w m (w m — 1), and thus, they 
move away from the point (uw, Z») = (1, 1), while they 
approach the point (w m *, Z*) = (1, — 1). In addition, the 
solution of this equation is u) m = [l + ce c( - 1+3,i ') r ]~ 1 , with 
c > an integration constant. Using this solution in 
equation H'/H = — Z(q+1) we find that for w = 1/3, 
H/H a — y/Lu m /(l — LJ m ), where H Q is another integra- 
tion constant. Then, the equation for iv m (t) becomes 
doj m /dt = —2tiLi m \J H^ui m — /i(l — ui m ) 2 , and can be inte- 
grated giving t as a function of uj m or H. Therefore, 
in the region of the big bang/big crunch singularity one 
obtains a(t) ~ ^/ei, p(t) ~£ -2 , as in the standard radia- 
tion dominated big-bang scenario. This means that for 
v < 3/2 the energy exchange has no observable effects 
close to the big bang/big crunch singularity. 




FIG. 1: Influx, w — 0, v < 3/2. The arrows show the direc- 
tion of increasing cosmic time. The dotted line corresponds 
to wde ~ —1. The region inside (outside) the dashed line 
corresponds to acceleration (deceleration). The region with 
Z > represents expansion, while Z < represents collapse. 
The present universe is supposed to be close to the global 
attractor. 

Since our proposal relies on the existence of an attrac- 
tor, we shall restrict ourselves to the case v < 3/2. It is 
convenient to discuss the four possible cases separately: 
(i) w — with influx. The generic behavior of the so- 
lutions of equations H50 M 51 fl ^ s shown in Figure 1. We 
see that all the expanding solutions approach the global 
attractor. Furthermore, there is a class of collapsing so- 
lutions which bounce to expanding ones. Finally, there 
are solutions which collapse all during their lifetime to 
a state with finite p and H. The physically interest- 
ing solutions are those in the upper part of the diagram 
emanating from the big bang (u, Z)fti(l, 1). These solu- 
tions start with a period of deceleration. The subsequent 
evolution depends on the value of 3(3 2 /\p\, which deter- 
mines the relative position of the dashed and dotted lines. 
Specifically, for 3(i 2 /\p,\ > 1 (the case of Figure 1) one dis- 
tinguishes two possible classes of universe evolution. In 
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FIG. 2: Outflow, w = 1/3, v < 3/2. The arrows show the 
direction of increasing cosmic time. The region inside (out- 
side) the dashed line corresponds to acceleration (decelera- 
tion). The region with Z > represents expansion, while 
Z<0 represents collapse. 



the first, the universe crosses the dashed line entering the 
acceleration era still with wde > — 1 , and finally it crosses 
the dotted line to wde < — 1 approaching the attractor. 
In the second, while in the deceleration era, it first crosses 
the dotted line to wde < — 1, and then the dashed line 
entering the eternally accelerating era. For 3/3 2 /|^t| < 1, 
the dotted line lies above the dashed line, and, conse- 
quently, only the second class of trajectories exists. To 
connect with the discussion in the introduction, notice 
that the Friedmann equation l|38|l can be written in the 
form (JIJ with dark energy pde = (fiip+fi) Using (1531) , 
the equation for puE takes the form (|2J with 



-1 



WDE' 



3(l-w m ) 



2Z 2 -w m -l-6(l-3u;) 



(3 2 u m (l-Z 2 ) 



Z 2 ~uj rl 



(55) 



The global attractor l|41|) - (|43|l satisfies relations 10 and 
consequently, wde evolves to the value Wde* given by 
Q. As for the bouncing solutions, they approach the 
attractor after they cross the line Z 2 — uj m , where wde 



jumps from +oo to — oo; however, the evolution of the 
observable quantities is regular. 

(ii) w = with outflow. The generic behavior in this case 
is obtained from Figure 1 by the substitution Z — > — Z 
and t — + — r, which reflects the diagram with respect to 
the uj m axis and converts attractors to repelers. 

(iii) w = l/3 with outflow. Figure 2 depicts the flow dia- 
gram of this case. Even though in the case of radiation in 
general wde > — 1/3 from equation J55j), there are both 
acceleration and deceleration regions. Furthermore, from 
equation 0} it is f2 m * > 1. 

(iv) ui=l/3 with influx. This arises like in (ii) by reflec- 
tion of Figure 2 and resembles Figure 1 . 



IV. OPEN QUESTIONS 

There is a number of features of the present scenario, 
which require further analysis, (a) One is the question 
of stability of the negative tension brane. In a fixed flat 
background, such a brane would be obviously unstable, 
by the formation of wild ripples on the brane. Such wild 
fluctuations increase the area of the brane and, with neg- 
ative tension, have energy unbounded from below. How- 
ever, here the situation is different, since the brane is in 
an AdS background and the effective cosmological con- 
stant on it is zero. The stability analysis of our scenario is 
currently under study, (b) It would be interesting to in- 
vestigate if the partial success of the present scenario per- 
sists after one tries to fit the supernova data, the detailed 
CMB spectrum and nucleosynthesis, (c) Finally, the 
construction of the complete higher dimensional theory, 
the specification of the nature of the content of the bulk 
and of the mechanism of energy exchange with the brane 
is another set of crucial open questions, which we hope 
to deal with in the not too distant future. 
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